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Abstract 

Modern ether-drift experiments in vacuum could in principle detect the tiny refractive in- 
dex that, in a flat-space picture of gravity, is appropriate for an apparatus placed on the 
Earth's surface. In this picture, in fact, if there were a preferred reference frame, light on the 
Earth would exhibit a slight anisotropy with definite quantitative differences from General 
Relativity. By re-analyzing the data published by two modern experiments with rotating 
optical resonators, and concentrating on the part of the signal that should be free of spurious 
systematic effects, we have found evidences that would support the flat-space scenario. 



1. Introduction 



The present generation of ether-drift experiments, using rotating optical resonators, is cur- 
rently pushing the relative accuracy of the measured frequency shifts to the level 0(1O -16 ). 
As we shall try to illustrate, this level of accuracy is crucial to obtain fundamental infor- 
mations on the space-time structure of the vacuum and the possible existence of a preferred 
reference frame. 

In this paper, we'll present a re-analysis of the observations reported in Ref. [I] for the 
anisotropy of the speed of light in the vacuum. This re-analysis, that improves with a higher 
statistics on our previous work [2] based on the data of Refs.[3l[I], leads to two conclusions: 
i) the experiments exhibit a non-zero daily average for the amplitude of the signal ii) the 
magnitude of this average amplitude is entirely consistent with the theoretical anisotropy 
parameter [SJ [SJ [7] 

l-Bthl ~ 3(AA vacuum - 1) ~ 42 ■ 1(T 10 (1) 

that, in the presence of a preferred frame, enters the two-way speed of light in the vacuum 
to 0{v 2 /c 2 ) 

mzM^B^t (2) 
c <r 

In Eq.([T]) A/"vacuum indicates the effective vacuum refractive index that is appropriate, in a 
flat-space picture of gravity, for an apparatus placed on the Earth's surface. 

We emphasize that the prediction Eq.([T]) refers to experiments performed in the vacuum, 
i.e. in the highest vacua attainable with present technology. Therefore, in principle, it should 
not be compared with other type of experiments as those of Ref. [8] where about 100% of the 
electromagnetic energy propagates within a dielectric medium with J\f ~ 3 [9] . 

After this general Introduction, the plan of the paper is as follows. In Sect. 2, we shall first 
illustrate the basic formalism and report the experimental data of Ref. pQ. Then, in Sect. 3, 
we shall use these published data to deduce a basic quantity of any ether-drift experiment: 
the daily average amplitude of the signal Aq that, so far, has never been reported by the 
experimental groups. Further, in Sect. 4, we shall try to estimate the possible systematic 
uncertainty on our value of Aq, by also comparing with Ref. [4]. Then, in Sect. 5, we shall point 
out that the experimental values of Aq for the two experiments, besides being well consistent 
with each other, are also in good agreement with the theoretical prediction expected in a 
flat-space picture of gravity, in the presence of a preferred reference frame. Finally, in Sect. 6, 
we shall present our summary and conclusions. 
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2. Basic formalism and experimental data 



The experimental data reported in Ref.pQ refer to 27 short-period observations, performed 
by the Berlin group [3] during a total period of 392 days. The starting point for our analysis 
is the expression for the relative frequency shift of two optical resonators at a given time t. 
For the Berlin experiment [3], this can be expressed as 

= S(t) sin 2uj Iot t + C(t) cos 2uj Iot t (3) 

where iv ro t is the rotation frequency of one resonator with respect to the other which is kept 
fixed in the laboratory and oriented north-south. We observe that the sine component of 
the signal was denoted as S(t) in Ref.[3] and as B(t) in Ref.pQ. Here, we shall maintain the 
notation of Ref.[3] where the Fourier expansions of S(t) and C(t) are expressed as 

S(t) = So + S s i sinr + S& cos r + S S 2 sin(2-r) + S C 2 cos(2t) (4) 

C(t) = C + C s % sin r + C c i cos t + C s2 sin(2r) + C c2 cos(2t) (5) 

r = u; s idt being the sidereal time of the observation in degrees and u> s id ~ 23^56' ' I n t r oducing 
the colatitude of the laboratory and the unknown average velocity, right ascension and 
declination of the cosmic motion with respect to a hypothetical preferred frame (respectively 
V, a and 7), one finds the expressions reported in Table I of Ref. [3] in the RMS formalism 

Co = -^(3cos2 7 -l) (6) 

1 1 

C s \ = -K sin 27 sin a sin 2% C c \ = -i^ sin 27 cos a sin 2^ (7) 

C S 2 = —K cos 2 7 sin 2a (1 + cos 2 x) = —K cos 2 7 cos 2a(l + cos 2 x) (8) 

V 2 

K = \B\^ (9) 



where 



is proportional to the light anisotropy parameter. The corresponding S— quantities are also 
given by (S = 0) 

Ssi = — ^- 5 cl = (10) 
cos x cos x 

_ 2cos X _ 2 cos x r ni x 

1 + cos^ x 1 + cos X 
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To compare with the similar Diisseldorf experiment of Ref.[4], one should just re-nominate 
the two sets 

(Co, C a i, C c i, C S 2, C C 2) — > (Co, C\, C2, C3, C4) (12) 

(So,S s i,Sci,S S 2,S C 2) — ► (Bq, Bi, B 2 , B 3 , B4) (13) 

and introduce an overall factor of two for the frequency shift since, in this case, two orthogonal 
cavities are maintained in a state of active rotation. 

As suggested by the same authors of Refs.pJl H], it is safer to concentrate on the ob- 
served time modulation of the signal, i.e. on the quantities C s i, C cl , C S 2, C C 2 and on their 
S-counterparts. In fact, the constant components Co and So = Bo are likely affected by spu- 
rious systematic effects. The experimental C& and Sk coefficients, as extracted from Fig. 2 
of Ref.pQ, are reported in our Tables 1 and 2. The quoted errors are both statistical and sys- 
tematical. For the quantities we are considering, according to Ref.[3|, this latter component 
should be very small. 



3. The daily average amplitude of the signal 

For our analysis, we shall re- write Eq.([3]) as follows 

= Alt) cos(2u rot t - 20 o (t)) (14) 

vo 

with 

C(t) = Ait) cos 20 o (t) Sit) = A(t) sin 29 (t) (15) 

9o (t) representing the instantaneous direction of a hypothetical ether-drift effect in the plane 
of the interferometer. 

Eq.([I3D, while fully equivalent to Eq.([3|), introduces in the analysis a positive-definite 
quantity, the amplitude of the signal A(t). The interest in a computation of A(t) can be easily 
understood by comparing the two elementary amplitudes C(t) and S(t) with the cartesian 
coordinates x and y of a 2-dimensional particle motion and A(t) with the radial coordinate 
r = y 1 x 1 + y 2 , so that x = r cos 9 and y = r sin 9. For definiteness let us denote by X{ ± Axj 
and yi±Ayi the individual measurements of the particle position and assume that, to a good 
accuracy, the average values are x = y = 0. It goes without saying that this situation does not 
determine all properties of the motion since there might be a plenty of non-trivial rotationally 
invariant motions that differ for the average radius. To determine this other parameter, one 
has to compute the radial distances = \ x 2 + y 2 for the individual measurements and 
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finally take their average. Of course, since now one deals with the positive-definite quantities 
7"j, their mean value f will be definitely different from \fx 2 + y 2 = and will depend on the 
experimental accuracy of the individual values rj db Arj. Assuming that the errors for the 
elementary entries Xj and y^ belong to a normal distribution, the errors Ar^ can be computed 
according to standard error propagation for composite observables, namely 

An = ^cos 2 O^Axi) 2 + sin 2 ei(A yi ) 2 (16) 

Therefore, the final answer about a possible non-zero f will depend on the confidence level 
of the point r = 0, i.e. on the distributions of the ratios Ri = rj/Arj. At the same time, 
when using the vectorial observables S(t) and C(t), noise tends to average out to zero for 
large sample of data. With a positive-definite quantity as A(t) this does not occur. Therefore 
one has to estimate the effect of noise that might mimic a true signal. This point will be 
discussed in Sect. 4 . 

After these preliminary considerations, let us return to our basic amplitude A(t) by ob- 
serving that, in the framework where the amplitudes S(t) and C(t) are represented as in 
Sect. 2, it can be expressed in terms of v(t), the magnitude of the projection of the cosmic 
Earth's velocity in the plane of the interferometer, and of the light anisotropy parameter as 

2 cr 

To compute v(t), we shall use the theoretical relations given by Nassau and Morse [llj 
from which A(t), S(t) and C{t) can be obtained up to an overall proportionality constant. 
These expressions are valid for series of short-period observations, as those performed in 
Refs.[H O Sj, where, within each experimental session, the kinematical parameters of the 
cosmic velocity V are not appreciably modified by changes in the Earth's orbital motion 
around the Sun. 

In this approximation, by introducing the magnitude V of the full Earth's velocity with 
respect to a hypothetic preferred frame S, its right ascension a and angular declination 7, 
we get 



cos z{t) = sin 7 sin <f> + cos 7 cos <f> cos(r — a) (18) 

sin z(t) cos 9o(t) = sin 7 cos 4> — cos 7 sin cj) cos(r — a) (19) 

sinz(t) sin#o(0 = cos 7 sin(r — a) (20) 

v(t) = Vsmz(t), (21) 
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where z = z(t) is the zenithal distance of V, <j) is the latitude of the laboratory and again 
r = u> s idt denotes the sidereal time of the observation in degrees. As one can check, by using 
the above relations, together with Eqs. (fT5|) and (fT7]) . and replacing \ = 90° — cj), one re-obtains 
the expansions for C(t) and S(t) reported in Eqs. ([!])- jSJ). 

After having obtained this first consistency check, we can now replace Eq. (|2ip into 
Eq. (JTTj) and, by adopting a notation of the type in Eqs.dH)-©, express A(t) as 



A(t) = A + A 1 sin r + A 2 cos r + A 3 sin(2r) + A 4 cos(2r) (22) 



where 



Aq = g 1-^1 ~~ sm2 7 cos2 X ~~ 2 cos2 7 s ^ n2 -0 

j4i = — — | if | sin 27 sin a sin 2^ A2 = — — \K | sin 27 cos a sin 2^ (24) 

1 1 
A3 = —-|if I cos 2 7 sin 2a sin 2 x = — -r\K | cos 2 7 cos 2a sin 2 x (25) 

On the basis of the above relations, it is now possible to extract the average amplitude Aq 
from the published data of Ref.pj. On the other hand, no information on the phase of the 
signal 9o(t) can be obtained. To this end, in fact, one would need the full amplitudes S(t) 
and C(t) at the sidereal times of the observations and these are not available. 

To obtain Aq, we observe that the daily averaging of the signal (here denoted by (..)), 
when used in Eq. (f22"]) produces the relation 

(A 2 (t)) = A 2 + ±(A 2 + A 2 + A 2 + A 2 ) (26) 

On the other hand, using Eqs.([4|), ([5]) and (fl~5l) . one also obtains 

( A 2 (t) ) = < C\t) + S 2 (t) ) = eg + s 2 + ±(C 2 U + s 2 u + c 2 22 + S 2 22 ) (27) 

where we have introduced the combinations 



Cn^JC 2 1+ C 2 cl c 22 = JC 2 2 + C 2 2 (28) 



S 11 ^JS 2 sl + S 2 1 s 22 ^JS 2 2 + S 2 c2 (29) 



As one can check, replacing the expressions (|23|) - ([25|) . Eq. (|26j) gives exactly the same result 
that one would obtain replacing the values for the C- and S- coefficients in Eq. (|27p . Therefore, 
one can combine the two relations and get 

A 2 (l + r) = C 2 + S 2 + \{C 2 n + s 2 u + e 2 2 + s 2 22 ) = e 2 + S 2 + Q 2 (30) 
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where 



and 



Q \ll(C 2 u + Sf 1 + Ci 2 + Sl 2 ) (31) 



' (Af + Al + Aj + Al) (32) 



2AI 



To evaluate Aq we shall proceed as follows. On the one hand, we shall compute the ratio 
r = r(7, x) using the theoretical expressions Eqs. (123[) - ()25p . This gives 

< r < 0.4 (33) 

for the latitude of the laboratories in Berlin [3] and Diisseldorf [4J in the full range < I7I < 
7r/2. On the other hand, we shall adopt the point of view of the authors of Refs.[TJ [3j 2] 
that, even when large non-zero values of Co and So are obtained (compare with the value 
Co = (—59.0 ± 3.4 db 3.0) • 10~ 16 of Ref.[4] and with the large scatter of the data reported in 
in Fig. 2 of Ref.[lJ or in Fig. 3 of Ref.[3]), tend to consider these individual determinations as 
spurious effects. 

This means to set in Eq. (|30p 

S = ( A(t)sm26 (t) ) ~0 (34) 

C = ( A(t) cos 20 o (i) ) ~ (35) 

We can thus define a daily average amplitude, say Aq, which is determined in terms of Cn, 
Sn, C22 and S22 alone as 

Aq = . = ~ (0.92 ± 0.08)Q (36) 
yl + r 

Here the uncertainty takes into account the numerical range of r in Eq. (|33|) . This value 
provides, in any case, a lower bound to its true experimental value since 



The data for the various coefficients are reported in our Table 3 together with the quantity 
Q. The errors on the various coefficients have been computed, according to standard error 
propagation for composite observables, starting from the errors for the basic entries C& and 
Sk reported in our Tables 1 and 2. 
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4. The value of Ao and its experimental uncertainty 

In this section we shall now try to deduce the average value of Aq (or more precisely its lower 
bound) for the Berlin experiment [U [3]. To this end, a first simple choice could be to take 
the weighted average Q of the 27 values of Q reported in Table 3 and then use Eq. ([36|) . This 
straightforward strategy gives 

Q = (13.0 ±0.7) • 1(T 16 (38) 

In the above relation the quoted error should be understood as purely statistical. In fact, 
by taking the weighted average of a large number N of measurements, one obtains an error on 
the mean that vanishes as 1/yN in the limit N — > oo. To estimate a possible systematic error 
to be added in Eq. (|38p we shall make the reasonable assumption that systematic uncertainties 
affect the value Qi of the i-th experimental session by an amount Q^ ys which is comparable 
to AQi. To get a numerical value, we shall follow three different methods that turn out to 
give similar results. 

In a first approach, one could estimate the systematic error on the mean by simply 
averaging the 27 determinations Q s ? s ~ AQi reported in Table 4. This gives the value 

~ Yi ? Qr ~ 44 ' 10 ~ 16 (39) 

i 

As a second approach, one can estimate a mean systematic component by using the weighted 
average expression 

In this case, by using again the data reported in Table 4, we obtain 

AQ {2) ~ 3.5 • 10" 16 (41) 

Finally, one can also consider the values of the signal and estimate the noise for the i-th 
experimental session from the ratio Ri = Qi/AQi. Taking the arithmetic mean of the 27 
determinations R{ reported in Table 4 gives a value R ~ 3.6 that, when using the weighted 
average Q = 13.0 • 10 -16 , gives the other estimate 

AQ( 3 ) „ £ ^ 3. 6 . i - 16 (42) 
R 

Therefore, by averaging the three different determinations in Eqs.(|39p. (|4ip and (|42p . one 
obtains the value 

AQ~3.8-10~ 16 (43) 
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that we shall take as our final estimate of the systematic error on the mean thus replacing 
Eq.flMD with 

Q = (13.0 ±0.7 ±3.8) • 1(T 16 (44) 
Prom this, by using Eq. ([36|) . we get 

Aq = (12.0 ± 1.0 ±3.5) • 10~ 16 Ref.p] (45) 

where the former error is due to the variation < r < 0.4 and the latter reflects our estimate 
of the noise. 

As a further control of the validity of our Eq . (|44|) . we report in Tables 5-7 the C and 
S coefficients of Ref.[l]. As one can check the value Q = (13.1 ± 2.1) • 10~ 16 is in excellent 
agreement with our result (|44p and corresponds to 

i = (12.1 ± 1.0 ± 2.1) -10~ 16 Ref.0] (46) 

in perfect agreement with (|45|) . 

5. An effective refractive index for the vacuum 

In this section, we shall point out that the two experimental values in Eqs. (|45p and (|46|) are 
well consistent with the theoretical prediction 

< = \\B^\°4 ~ (9.7 ±3.5) • 10- 16 (47) 
2 c z 

of Refs.[6l[7]. This was obtained, in connection with the anisotropy parameter [5] |-E?th| ~ 
42 • 10 -10 , after inserting the average cosmic velocity (projected in the plane of the interfer- 
ometer) vq = (204 ±36) km/s that derives from a re-analysis [6j[7] of the classical ether-drift 
experiments. Due to this rather large theoretical uncertainty, the different locations of the 
various laboratories and any other kinematical property of the cosmic motion can be neglected 
in a first approximation. 

For a proper comparison, we also remind that in Refs.jH [7] the frequency shift was 
parameterized as 

^l = \B th \ V 4co S 26 (48) 

This relation is appropriate for a symmetrical apparatus with two rotating orthogonal lasers, 
where the difference c{9) — c(0) gets replaced by c(tt/2 + 9) — c(9) as in the Diisseldorf 
experiment [3], and gives an average amplitude 

A s ymm = 2A t h ~ (19 ± 7) • 10- 16 (49) 
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The theoretical prediction for the anisotropy parameter was obtained starting from the 
many analogies that one can establish, in the weak-field limit, between General Relativity 
and a flat-space description where gravity leads to re-defined masses, space-time units and 
an effective vacuum refractive index. This alternative approach, see for instance Wilson |12j . 
Gordon [13], Rosen p3], Dicke [15] . Puthoff [16] and even Einstein himself |17j . before his 
formulation of a metric theory of gravity, in spite of the deep conceptual differences, produces 
an equivalent description of the phenomena in a weak gravitational field. 

The substantial phenomenological equivalence of the two approaches was well summarized 
by Atkinson as follows [18] : "It is possible, on the one hand, to postulate that the velocity of 
light is a universal constant, to define natural clocks and measuring rods as the standards by 
which space and time are to be judged and then to discover from measurement that space- 
time is really non-Euclidean. Alternatively, one can define space as Euclidean and time as 
the same everywhere, and discover (from exactly the same measurements) how the velocity 
of light and natural clocks, rods and particle inertias really behave in the neighborhood of 
large masses." 

This type of analogy, which is preserved by the weak-field classical tests, is interesting in 
itself and deserves to be explored. In fact, "...it is not unreasonable to wonder whether it 
may not be better to give up the geometric approach to gravitation for the sake of obtaining 
a more uniform treatment for all the various fields of force that are found in nature" |14j . 

To distinguish between the two interpretations, one can start from the zeroth order ap- 
proximation to the problem embodied in the Equivalence Principle. As it is well known, 
this basic principle, introduced before General Relativity, does not necessarily rely on the 
notion of a curved space-time (see e.g. Ref.[19] or even the critical point of view of Ref.[20j). 
According to it, for an observer placed in a freely falling frame, local Lorentz invariance is 
valid to first order. Therefore, given two space-time events that differ by (dx , dy , dz , dt) , and 
the space-time metric 



one gets from ds 2 = the same speed of light that one would get in the absence of any 
gravitational effect. 

However, to a closer look, an observer placed on the Earth's surface is equivalent to a 
freely-falling frame up to the presence of the Earth's gravitational field. In this case, both 
General Relativity and the flat-space approach predict the weak-field, isotropic form of the 




(50) 
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metric 



ds 2 = c 2 dt 2 (l - ^) - (1 + ^^)(dx 2 + dy 2 + dz 2 ) = c 2 dr 2 - dl 2 (51) 
c A R c z R 



G being Newton's constant and M and R the Earth's mass and radius. Here dr and dl denote 
respectively the elements of "proper" time and "proper" length in terms of which, in General 
Relativity, one would again deduce from ds 2 = the same universal value 4- = c. 

On the other hand, in the flat-space approach the condition ds 2 = is interpreted in 
terms of an effective refractive index for the vacuum 

2GM m , . 

AA V acuum - 1 ~ ~ 14 • 10- 10 (52) 

as if Euclidean space would be filled by a very rarefied " gravitational medium" (let us recall 
that a moving dielectric medium acts on light as an effective gravitational field, see for 
instance Refs.pSl EI]). 

Now, since in the flat-space approach light can be seen isotropic by only one inertial frame 
[22], say E, the ether-drift experiments can clarify whether £ coincides with the Earth's frame 
or with a hypothetical preferred frame. In the former case, corresponding to no anisotropy 
of the two-way speed of light in the vacuum, the equivalence between General Relativity and 
the gravitational-medium picture would persist. In the latter case, for the observer sitting 
on the Earth, there would be off-diagonal elements in the metric 

90i ~ (A/" v 2 acuum - 1) (53) 

and an anisotropy of the speed of light. Its value can be estimated directly in flat space, by 
Lorentz transforming from the isotropic value c/A/" va cuum m S, and one predicts an anisotropy 
of the two-way speed of light 

m=m~Bj***, (54, 

c c z 

governed by the parameter El [7] 

|Sth|~3(AA„ m -l)~ 42-10- 10 (55) 

For this reason, assuming a typical value V 2 /c 2 ~ 10 -6 , the present ether-drift experiments, 
with their 0(1O -16 ) accuracy, represent precision probes of the vacuum and of its space-time 
structure. In this respect, it is interesting that, independently of our re-analysis of the data, 
a 2 — 3 a non-zero signal in the photon sector h _ = (—10.3 ± 3.9) • 10~ 16 (see Table II of 
Ref.[l]) was reported by the authors of Ref.p] by using the SME parametrization [23]. 
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6. Summary and conclusions 



In this paper, we have presented a re-analysis of the Berlin [HE] and Diisseldorf [2] ether-drift 
experiments that, by employing rotating optical resonators, attempt to establish the isotropy 
of the speed of light in the vacuum to a level of accuracy 0(1O -16 ). In our re-analysis we have 
extracted from the published data a basic observable of an ether-drift experiment, namely 
the daily average for the amplitude of the signal. By denoting with (..) the daily average of 
any given quantity, this is defined as Aq = (A(t)) after re-writing the frequency shift 

Au(t) 



in the form 



Sit) sin 2uj Iot t + C (t) cos 2w ro ti (56) 



= A(t) cos(2u rot t - 20 o (t)) (57) 

By further assuming, as the authors of Refs.[3l|4j do, that experimental results providing large 
non-zero values for either (C(t)) = Co or {S(t)) = So in Eqs.flU) and ([5]) should be interpreted 
as spurious effects (e.g. due to thermal drift, non-uniformity of the rotating cavity speed, 
misalignment of the cavity rotation axis,...), Aq can be approximated by its lower bound 

Aq ~ (0.92 ± 0.08)Q (58) 

where 



Q" \/^^i 2 i + 5 1 2 1 + C| 2 + S 2 2 2 ) (59) 



is given in terms of the coefficients Cn, C22, Sii, S22 defined in Eqs. (j28]) - (|29j) . These coef- 
ficients, that reflect the time modulation of the signal, are much less affected by spurious 
effects (as compared to Co and So) and so should be the value of Aq. 

Now, the two resulting experimental determinations in Eqs. (|45p and (06]), namely Aq = 
(12.0 ± 1.0 ± 3.5) • 10~ 16 from Ref.p] and i = (12.1 ± 1.0 ± 2.1) • 10~ 16 from Ref.0], 
besides being in excellent agreement with each other (and with the previous lower-statistics 
re-analysis of the Berlin experiment reported in [2j) also agree with the theoretical prediction 
(9.7 ± 3.5) • 10~ 16 of Refe.0 E]- This prediction was obtained in a flat-space description of 
gravity, by assuming an average speed with respect to a hypothetical preferred frame of the 
type suggested by a re-analysis of the classical ether-drift experiments. Therefore this non- 
trivial level of consistency, between different experiments and with a theoretical prediction 
formulated before the experiments were performed, supports the conclusion that a non-zero 
anisotropy of the speed of light in the vacuum has actually been measured with values of 
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the anisotropy parameter that are one order of magnitude larger than those quoted by the 
authors of Refs.^IJ. 

A first obvious possibility for this discrepancy is that in our average we might have grossly 
underestimated the effect of the residual noise, even in the time- dependent part of the signal. 
In this case, our error in Eq. (|44p would be too small. For instance, by insisting on a zero 
signal and deciding that a vanishing amplitude requires consistency with the value Q = 
to some definite confidence level, say 10% one side, one might argue that the correct final 
estimate for the Berlin experiment should be Q ~ (13 ± 10) • 10~ 16 . To this end, however, 
the individual determinations in Table 4 should exhibit a typical value Ri = Qi/AQi ~ 1.3 
which is quite different from their average result Ri ~ 3.6. Analogous considerations apply 
if the confidence level were pushed to 5%, again one side. In this case, the final estimate 
for the Berlin experiment should be Q ~ (13 db 8) • 10~ 16 with individual determinations 
Ri = Qi/AQi ~ 1.6 in Table 4. 

Now, the individual central values Qi and their errors AQi reported in Table 4 were 
computed, starting from the values of the elementary entries and Sk reported in Tables 1 
and 2, according to the standard statistical rules for composite observables. Errors such as 
those needed to obtain Q ~ (13 ± 10) • 10~ 16 or Q ~ (13 ±8) • 10~ 16 would be respectively five 
or four times larger than that extracted from the data of Ref . [4] after a single experimental 
session (see Table 7). This would be equivalent to deny the assumed high precision of the 
Berlin experiment. 

On the other hand, by accepting our estimate of the errors, the reasons for the discrepancy 
with the anisotropy parameter reported in Refs.[3l 0] have to be searched elsewhere. To 
this end, we observe preliminarily that introducing C\i, C22, S\\ and S22 represents the 
simplest way to obtain rotationally invariant combinations, out of the elementary coefficients 
C s i, C c i.. and of their S-counterparts. Their use can eliminate possible spurious effects, 
depending on the relative phases, that enter the delicate splitting of the signal in its various 
Fourier components (see the discussion given by the authors of Ref. [3] and their note [13]). 

Now, suppose that for some reason, during the i-th experimental session the coefficients 
C s i and C c i were erroneously rotated by an angle <pi with respect to their true value. This 
means erroneous values 



C s i(i) = C s i cos ^ - C c i sin ipi 



(60) 



and 



C c i{i) = C c i cos ifi + C s i sin tpi 



(61) 
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with analogous effects for the other parameter pairs (C S 2,C C 2), (S s i,S c i) and (S C 2,S S 2) and 
an overall shift in the value of the right ascension a — > a + ^ . 

Nevertheless, even if this happens, our rotationally invariant combinations would remain 
unchanged, i.e. 



with similar results for C22, Six and S22 ■ Therefore, even substantial levels of phase error, 
that can produce zero intersession averages for the various Ck and Sk vectorial coefficients 
(and thus a vanishing anisotropy parameter from these averages), would not change our 



apply to the separation of the signal in its S(t) and C{t) parts which is obtained by extracting 
the components proportional to sin 2uj TO tt and cos 2u! ro tt. Again, the combination A 2 = 
S 2 + C 2 is not modified by replacing io rot t — ► u; ro ti + fi$. This type of argument would also 
explain the remarkable agreement between our final estimate Eq . (|44[) from Ref.pQ and the 
value Q = (13.1 ± 2.1) • 10~ 16 from Ref.[3], in spite of the completely different Ck and Sk 
coefficients obtained by the two experiments. 

Clearly, the simplest way to check our result would be that the authors of Refs.[TJ [3j 2] 
could repeat their analysis of the data, replacing Eq.© with Eq.(fl4"]). and compute from 
scratch Aq (and its lower bound Aq) for the various experimental sessions. This computation, 
that would only require the elementary algebraic relations used in this paper, could also 
provide a powerful consistency check of the whole experiment. At the present, since this has 
not been done, by accepting the errors of Refs.[U [3j H] for the elementary entries C s i, C c \... 
and their S-counterparts, our values of Aq are the only existing estimate of this basic physical 
quantity. 

For this reason, since with our analysis one would be driven to deduce non-trivial conse- 
quences, such as the existence of a preferred frame and a flat-space description of gravity, we 
emphasize the importance of comparing different points of view and approaches to the data 
to finally achieve a full understanding of the underlying fundamental physical problem. 
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Table 1: The experimental Ck coefficients as extracted from Fig. 2 of Ref.[lJ. 



C,i[xl0- 16 ] 


C cl [xl0" 16 ] 


C s2 [xW- 16 } 


C c2 [xl0- 16 ] 


9.2 ± 4.6 


10.0 ± 5.4 


— 1.5 ± 4.6 


0.0 ± 5.4 


—2.3 ± 6.9 


—4.6 ± 6.2 


-5.8 ± 6.5 


-0.8 ± 6.2 


—2.3 ± 3.8 


0.0 ± 3.1 


—2.3 ± 3.8 


—4.6 ± 3.8 


4.6 ± 4.6 


11.5 ± 3.8 


-14.6 ± 3.8 


-0.8 ± 4.6 


5.4 ± 8.5 


-11.5 ± 8.5 


11.5 ± 8.5 


-4.6 ± 6.9 


—5.4 ± 5.4 


—6.2 ± 6.2 


6.9 ± 5.4 


1.5 ± 5.4 


— 10.8 ± 6.2 


—5.4 ± 5.4 


—6.2 ± 6.2 


—4.6 ± 5.4 


—3.1 ± 9.2 


7.7 ± 9.2 


— 10.8 ± 9.2 


7.7 ± 8.5 


-6.2 ± 6.2 


6.2 ± 7.7 


-2.7 ± 6.5 


6.9 ± 6.2 


— 1.5 ± 4.6 


—9.2 ± 4.6 


— 10.0 ± 5.4 


-3.1 ± 5.4 


2.3 ± 3.8 


4.6 ± 3.1 


-3.1 ± 3.1 


0.0 ± 3.8 


13.8 ± 7.7 


— 10.0 ± 6.9 


10.0 ± 6.9 


3.8 ± 7.7 


-20.0 ± 4.6 


-5.4 ± 5.4 


4.6 ± 4.6 


9.2 ± 5.4 


-2.3 ±3.8 


6.2 ±3.1 


-3.8 ±3.8 


1.5 ±3.8 


-16.2 ±6.9 


22.3 ±6.9 


-3.1 ±6.2 


10.0 ±6.2 


-6.2 ±7.7 


18.5 ±9.2 


6.2 ±9.2 


13.1 ±7.7 


3.1 ±3.1 


3.8 ±3.8 


-10.8 ±3.1 


5.4 ±3.1 


-10.0 ±8.5 


1.5 ±7.7 


-1.5 ±7.7 


3.1 ±8.5 


0.8 ±8.5 


6.9 ±8.5 


-13.1 ±8.5 


2.3 ±9.2 


-16.2 ±6.9 


4.6 ±7.7 


0.8 ±6.9 


26.2 ±8.5 


1.5 ±6.2 


-21.5 ±6.2 


12.3 ±6.2 


-6.2 ±5.4 


0.8 ±6.9 


-1.5 ±7.7 


-7.7 ±7.7 


-13.1 ±7.7 


7.7 ±9.2 


-14.6 ±8.5 


14.6 ±8.5 


3.8 ±9.2 


0.0 ±4.6 


6.9 ±3.8 


2.3 ±3.8 


7.7 ±3.8 


-16.2 ±6.9 


-10.0 ±8.5 


1.5 ±7.7 


11.5 ±7.7 


-1.5 ±4.6 


0.8 ±3.8 


-3.1 ±4.6 


1.5 ±3.8 


-5.4 ± 10.0 


10.8 ± 12.3 


-20.8 ± 10.0 


-14.6 ± 10.8 
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Table 2: The experimental Sk = Bk coefficients as extracted from Fig. 2 of Ref.JT]. 



^i[xl0" 16 ] 


S c i[xl0" 16 ] 


MxlO" 16 ] 


5 c2 [xl0" 16 ] 


2.3 ± 3.8 


— 10.8 ± 4.6 


—4.6 ± 4.6 


—4.6 ± 4.6 


14.6 ± 6.9 


-10.0 ± 6.9 


—6.9 ± 6.9 


1.5 ± 6.2 


— 1.5 db 4.6 


—5.8 ± 3.5 


3.8 ± 3.8 


—2.3 ± 3.8 


-8.5 ± 3.8 


8.5 ± 3.8 


-10.0 ± 3.8 


4.6 ± 4.6 


3.1 ± 7.7 


-3.8 ± 8.5 


-0.8 ± 6.9 


-6.9 ± 6.9 


1.5 ± 6.2 


3.1 ± 6.2 


— 10.8 ± 6.2 


1.5 ± 6.2 


9.2 ± 6.2 


3.8 ± 5.4 


5.4 ± 5.4 


— 1.5 ± 6.2 


—3.8 db 8.5 


—2.3 ± 8.5 


— 16.2 ± 8.5 


— 14.6 ± 8.5 


—4.6 d= 6.2 


10.0 ± 6.9 


9.2 ± 6.2 


1.5 ± 6.2 


3.8 ± 5.4 


—4.6 ± 4.6 


5.8 ± 5.0 


4.6 ± 4.6 


3.1 ± 3.1 


1.5 ± 3.1 


-3.8 ± 3.8 


-1.5 ± 3.1 


3.1 ± 7.7 


2.3 ± 6.9 


2.3 ± 6.9 


0.8 ± 6.9 


-3.8 d= 3.8 


0.4 ± 5.0 


3.8 ± 3.8 


-2.3 ± 3.8 


5.4 ±3.8 


1.5 ±3.1 


6.2 ±3.1 


3.1 ±3.1 


-14.6 ± 5.4 


20.0 ±6.2 


3.8 ±5.4 


6.5 ±5.8 


10.0 ±6.9 


-19.2 ±8.5 


16.9 ±7.7 


-16.2 ±6.9 


-2.3 ±3.8 


2.3 ±3.8 


8.5 ±3.8 


-8.5 ±3.8 


2.3 ±6.9 


-10.0 ±6.9 


-8.5 ±6.9 


1.5 ±7.7 


-10.0 ±8.5 


3.8 ±8.5 


-5.4 ±8.5 


-3.1 ±9.2 


0.0 ±7.7 


6.2 ±7.7 


10.8 ±7.7 


-4.6 ±9.2 


2.3 ±5.4 


-29.2 ±6.2 


9.2 ±6.2 


-8.5 ±5.4 


-6.9 ±6.9 


-0.8 ±6.9 


-11.5 ±6.9 


-16.2 ±6.9 


-6.2 ±9.2 


-7.7 ±9.2 


10.8 ±9.2 


-4.6 ±9.2 


5.4 ±3.8 


8.5 ±3.8 


-8.5 ±3.8 


0.8 ±3.8 


7.7 ±7.7 


-1.5 ±9.2 


-8.5 ±8.5 


5.4 ±8.5 


-0.8 ±3.8 


-0.8 ±3.8 


-2.3 ±3.8 


3.8 ±3.8 


-16.2 ± 11.5 


6.2 ± 13.8 


-25.4 ± 11.5 


-16.2 ± 11.5 
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Table 3: The experimental values of Ref.Jlj for the combinations of C— and S— coeffi- 
cients defined in Eqs. (l28p - (|29p and the resulting Q from Eq. (|3ip , For simplicity, we report 
symmetrical errors. 



Cn[xl0- 16 ] 


C 22 [xl0- 16 ] 


MxlO- 16 ] 


<S 22 [xlO- 16 ] 


Q[xl0" 16 ] 


13.6 ± 5.0 


1.5 ± 4.6 


11.0 ± 4.6 


6.5 ± 4.6 


-it) t) 1 o A 

13.3 ± 3.4 


5.2 ± 6.3 


5.8 ± 6.5 


17.7 ± 6.9 


7.1 ± 6.9 


1 A C 1 A O 

14.6 ± 4.8 


O O 1 no 

2.3 ± 3.8 


5.2 ± 3.8 


6.0 ± 3.5 


4.5 ± 3.8 


6.6 ± 2.6 


1 ft i | A f\ 

12.4 ± 4.0 


14.6 ± 3.8 


12.0 ± 3.8 


11.0 ± 4.0 


"1 ^7 O 1 o o 

17.8 ± 2.8 


12.7 ± 8.5 


-t ft A | O O 

12.4 ± 8.3 


4.9 ± 8.2 


7.0 ± 6.9 


14.0 ± 5.8 


8.1 ± 5.9 


7.1 ± 5.4 


3.4 ± 6.2 


10.9 ± 6.2 


11.1 ± 4.2 


12.0 ± 6.0 


7.7 ± 5.9 


10.0 ± 6.0 


5.6 ± 5.4 


13.0 ± 4.2 


8.3 ± 9.2 


■i c% ft i r\ r\ 

13.2 ± 9.0 


4.5 ± 8.5 


21.8 ± 8.5 


19.2 ± 6.1 


8.7 ± 7.0 


7.4 ± 6.2 


11.0 ± 6.8 


9.4 ± 6.2 


-It) 1 \ 1 A *~7 

13.0 ± 4.7 


9.4 ± 4.6 


10.5 ± 5.4 


6.0 ± 4.9 


7.4 ± 4.9 


12.0 ± 3.5 


5.2 ± 3.2 


3.1 ± 3.1 


3.4 ± 3.1 


4.1 ± 3.7 


5.7 ± 2.4 


17.1 ± 7.4 


10.7 ±7.0 


3.8 ± 7.4 


2.4 ±6.9 


14.6 ± 5.2 


20.7 ±4.7 


10.3 ±5.2 


3.9 ±3.9 


4.5 ±3.8 


16.9 ±3.3 


6.6 ±3.2 


4.1 ±3.8 


5.6 ±3.8 


6.9 ±3.1 


8.3 ±2.4 


27.5 ± 6.9 


10.5 ±6.2 


24.8 ±5.9 


7.6 ±5.7 


27.7 ±4.5 


19.5 ±9.1 


14.5 ±8.0 


21.7 ±8.2 


23.4 ± 7.3 


28.3 ±5.7 


4.9 ±3.6 


12.0 ±3.1 


3.3 ±3.8 


12.0 ±3.8 


12.7 ±2.5 


10.1 ±8.4 


3.4 ±8.3 


10.3 ±6.9 


8.6 ±6.9 


12.1 ± 5.3 


7.0 ±8.5 


13.3 ±8.5 


10.7 ±8.5 


6.2 ±8.7 


13.7 ±6.0 


16.8 ±7.0 


26.2 ±8.5 


6.2 ±7.7 


11.7 ±8.0 


23.9 ±5.7 


21.6 ±6.2 


13.8 ±6.0 


29.3 ±6.1 


12.5 ±5.8 


28.9 ±4.3 


1.7 ±7.5 


15.2 ±7.7 


7.0 ±6.9 


19.9 ±6.9 


18.4 ± 5.1 


16.5 ±8.6 


15.1 ±8.5 


9.9 ±9.2 


11.7 ±9.2 


19.2 ±6.2 


6.9 ±3.8 


8.0 ±3.8 


10.0 ±3.8 


8.5 ±3.8 


11.9 ±2.7 


19.0 ± 7.4 


11.6 ±7.7 


7.8 ± 7.8 


10.0 ±8.5 


18.1 ± 5.4 


1.7 ±4.5 


3.4 ±4.5 


1.1 ±3.8 


4.5 ±3.8 


4.2 ±2.9 


12.0 ± 11.9 


25.4 ± 10.3 


17.3 ± 11.9 


30.1 ± 11.5 


31.6 ±7.9 
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Table 4: For each experimental session, we report the values Qi, AQi and the ratio Ri = 
Qi/AQi from Table 3. 



Qi[xW iD ] 


A /I r -i r\ — 1 fS 1 

AQj[xlO iD ] 


Ri = Qi/AQi 


13.3 


3.4 


3.9 


14.6 


4.8 


3.0 


6.6 


2.6 


2.5 


17.8 


2.8 


6.3 


14.0 


5.8 


2.5 


11.1 


4.2 


2.6 


13.0 


4.2 


3.1 


19.2 


6.1 


3.1 


13.0 


4.7 


2.8 


12.0 


3.5 


3.4 


5.7 


2.4 


2.4 


14.6 


5.2 


2.8 


16.9 


3.3 


5.1 


8.3 


2.4 


3.4 


27.7 


4.5 


6.2 


28.3 


5.7 


5.0 


12.7 


2.5 


5.1 


12.1 


5.3 


2.3 


13.7 


6.0 


2.3 


23.9 


5.7 


4.2 


28.9 


4.3 


6.7 


18.4 


5.1 


3.6 


19.2 


6.2 


3.1 


11.9 


2.7 


4.4 


18.1 


5.4 


3.3 


4.2 


2.9 


1.4 


31.6 


7.9 


4.0 
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Table 5: The experimental C— coefficients as reported in Ref.|4 



C,i[xl0- 16 ] C c i[xl0- 16 ] C s2 [xl0- 16 ] C c2 [xl0" 16 ] 
-3.0 ±2.0 11.0 ±2.5 1.0 ± 2.5 0.1 ± 2.5 
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Table 6: The S— coefficients of Ref.[4|. The values for the S-coefficients, constrained by the 
authors of Ref.[4] in their fits to the data to the theoretical predictions Egs- fTUj) and (fTTj) . 
have been deduced from Table 5 using these relations. 



S s i[xl0- 16 ] 5 c i[xl0- 16 ] 5 s2 [xl0- 16 ] 5 c2 [xl0" 16 ] 
-14.1 ±3.2 -3.8 ±2.6 -0.1 ±2.5 1.0 ±2.5 
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Table 7: The values of Ref.[l] for the combinations of C— and S— coefficients defined in 
Eqs. ([25]) - ([2U]) and the resulting Q defined in Eq. ([3"Tj) . as computed from Tables 5 and 6. For 
simplicity, we report symmetrical errors. 



CiibdO" 16 ] 


C 22 [xl0- 16 ] 


Sn[xl0- 16 ] 


5 22 [xl0- 16 ] Q[xl0- 16 ] 


11.4 ±2.5 


1.0 ±2.5 


14.6 ±3.3 


1.0 ±2.5 13.1 ±2.1 
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